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Abstract 

We study the implications of thermal leptogenesis for neutrino parameters. As- 
suming that decays of A'^i, the lightest of the heavy Majorana neutrinos, initiate 
baryogenesis, we show that the final baryon asymmetry is determined by only four 
parameters: the CP asymmetry ei, the heavy neutrino mass Mi, the effective light 
neutrino mass mi, and the quadratic mean m of the light neutrino masses. Impos- 
ing the CMB measurement of the baryon asymmetry as constraint on the neutrino 
parameters, we show, in a model independent way, that quasi-degenerate neutri- 
nos are incompatible with thermal leptogenesis. For maximal CP asymmetry £i, 
and neutrino masses in the range from (Am^^j)^/^ to (Am^^^)"^/^, the baryogenesis 
temperature is Tb = ©(lO^") GeV. 



1 Introduction 



The explanation of the cosmological baryon asymmetry is a challenge for particle physics. 
In principle, already the standard model contains all necessary ingredients, baryon number 
violation, C and CP violation, and also the required departure from thermal equilibrium 
could be generated during the electroweak phase transition [1]. However, due to the lower 
bound on the Higgs boson mass from LEP, electroweak baryogenesis is no longer a viable 
mechanism, except for some supersymmetric extensions of the standard model [2]. 

A simple and elegant explanation of the observed baryon asymmetry is offered by 
neutrino physics. During the past years data on atmospheric and solar neutrinos have 
provided strong evidence for neutrino masses and mixings. In the seesaw mechanism 
[3] the smallness of these neutrino masses is naturally explained by the mixing of the 
left-handed neutrinos with heavy Majorana neutrinos. Further, the connection between 
baryon and lepton number in the high-temperature, symmetric phase of the standard 
model due to rapid sphaleron transitions [4] is by now firmly established [5]. As in classical 
GUT baryogenesis [6] , out-of-equihbrium decays of the heavy Majorana neutrinos can then 
generate a lepton asymmetry which, by sphaleron processes, is partially transformed into 
a baryon asymmetry [7] . 

A beautiful aspect of this 'leptogenesis' mechanism is the connection between the 
cosmological baryon asymmetry and neutrino properties. This connection is established 
by standard kinetic calculations [8, 9], very much like in big bang nucleosynthesis [6], 
where light nuclei play the role analogous to leptons in leptogenesis. The requirement of 
'successful baryogenesis', i.e. the existence of neutrino masses and mixings for which the 
predicted and the observed value of the baryon asymmetry are in agreement, constitutes 
a severe test for models of neutrino masses, which has been extensively explored during 
the past years [10]. 

On the experimental side, the precision of measurements of the baryon asymmetry has 
significantly improved with the observation of the acoustic peaks in the cosmic microwave 
background radiation (CMB). The BOOMERanG and DASI experiments have measured 
the baryon asymmetry with a (Icr) standard error of ~ 15% [11, 12], 

(QB/^^)^^^ = 0.022in^^ (1) 

Since the number of relic photons per comoving volume is very unlikely to have changed 
after recombination, this is easily translated into a measurement of the quantity tjb = 
{nB/uj) at the present time, 

ri^sr - (e.OlSi) X lO-^-^ . (2) 
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In the near future the MAP experiment will provide results of the first full sub-degree 
sky survey of temperature anisotropies [13]. The expected (Icr) standard error on rjBo 
is ~ 10%, using only CMB data, no polarization measurement and even allowing for 
the presence of gravity wave perturbations [14]. The PLANCK satellite, whose launch 
is planned for 2007 [15], should reduce this error to ~ 1%. If the polarization will be 
measured, and if it will be possible to add extra CMB information on other cosmological 
parameters, then the error may become even less than 1% [14]. 

In the following we shall calculate the baryon asymmetry by solving the Boltzmann 
equations given in [8, 9], assuming that the dominant contribution is given by decays of 
Ni, the lightest of the heavy Majorana neutrinos. This assumption is well justified in the 
case of a mass hierarchy among the heavy neutrinos, i.e. Mi ^ M2. Mo,, and it is also 
known to be a good approximation, if M2.3 — Mi = 0{Mi) [16]. The case M2,3 — Mi <^ Mi 
requires a special treatment. For some flavour structures of the neutrino mass matrices it 
is also conceivable that the decays of the heavier neutrinos N2 or A^3 are the main source 
of the baryon asymmetry [17]. In our analysis wc shall assume that quantum corrections 
[18] to the Boltzmann equations are small. So far, no detailed quantitative study of this 
important question has been carried out. 

As we shall see, within this framework one is left with only four parameters: Mi, the 
CP asymmetry ei, the effective neutrino mass rhi and m, the quadratic mean of the light 
neutrino masses. For each set of values of these four parameters the Boltzmann equations 
yield a prediction for the baryon asymmetry. A comparison with the observed value then 
defines an allowed region in the space of neutrino parameters. 

In addition to the neutrino parameters, there are also three quantities which character- 
ize the initial conditions: the initial temperature, the initial abundance of heavy neutrinos 
and, of course, the initial baryon asymmetry. A detailed study of the stability of the final 
baryon asymmetry under variations of these initial conditions will be presented elsewhere 
[16]. In the following we shall illustrate this dependence by presenting all results for two 
different choices of the initial A^i abundance, namely zero and thermal initial abundance. 
All values in-between can be estimated by interpolation between these two cases. One 
can also easily extrapolate the results to initial A^i abundances higher than the thermal 
one. Fortunately, for the most interesting range of rhi the dependence on the initial Ni 
abundance turns out to be very small. Clearly, a theory of the very early universe, like 
inflation, is needed to calculate the initial conditions for baryogenesis. 

The paper is organized as follows. In section 2 we recall the Boltzmann equations which 
we then solve numerically. We also briefly discuss some approximations underlying these 
equations. Section 3 deals with the theoretically allowed range of the neutrino parameters. 
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in particular the upper bound on the CP asymmetry. In section 4 we present our numerical 
results for the baryon asymmetry and discuss the dependence on the neutrino parameters 
and on the choice of the initial condition. We then investigate the constraints imposed 
by the CMB result on the neutrino parameters. Our conclusions are given in section 5. 

2 Solutions of the Boltzmann equations 

The dynamical generation of a baryon asymmetry requires that the particle interactions 
do not conserve baryon number, C and CP. In leptogenesis these conditions are realized 
by the couplings of the heavy Majorana neutrinos Ni. Their decays can generate an 
asymmetry in the number of leptons and antileptons, and therefore in B — L. The 
crucial departure from thermal equilibrium is provided by the expansion of the universe. 
At temperatures T ~ 0{Mi) the abundance of heavy neutrinos exceeds the thermal 
abundance due to their weak interactions with the thermal bath. 

A quantitative description of this non-equilibrium process is obtained by means of 
kinetic equations. The relevant processes in the thermal plasma are: 

• A^i decays (D) and inverse-decays (ID) into leptons and Higgs bosons, A^i <-> (f)l, 
and into antileptons and anti-Higgs bosons, Ni ^l; 

• AL = 2 scatterings mediated by the exchange of all heavy Majorana neutrinos, 
14> ^ 14> (iV), and // ^ 4>4>, U^4>(j) {N, t); 

• AL = 1 scatterings, Nil(l) ^ t(t)q(q) (0, s) and Nit(t} ^ l(l)q(q) (0,^); 

in brackets we have indicated how the rates of these processes are labeled in the following. 
In principle, one could also have additional processes, in particular those which contribute 
to bring the heavy neutrinos initially into thermal equilibrium [9]. In the present, minimal 
framework we neglect such interactions. 

Since we are assuming that A^i decays are the origin of lepton and baryon asymmetries, 
the natural temperature scale is given by the mass Mi. It is therefore convenient to 
measure temperature in units of Mi and to introduce the dimensionless variable z — 
Mi/T. For realistic values Mi ^ 100 GeV, all standard model particles can be treated as 
massless, and we shall assume that they are in thermal equilibrium. 

The time evolution of a charge density or a number density nx depends on the micro- 
physical processes in the thermal plasma as well as the expansion of the universe. For the 
discussion of leptogenesis it is convenient to consider instead of the number density nx 
the particle number Nx in some portion of comoving volume, which takes the effect of 
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the expansion automatically into account. We choose the comoving volume which 
contains one photon at time t^, before the onset of leptogenesis, 



Nx{t)=nx{t)R,{tf , 



(3) 



with 



-1/3 



(4) 



and therefore A^-y(t^) = 1. For a boson with qb degrees of freedom one has A''^(t^) = Qb/'^i 
whereas for a fermion with g^,- degrees of freedom Np(t^) = Sgp/S. Alternatively, one 
may normalize the number density to the entropy density s and consider Yx = nx/s, as 
frequently done in the literature. If entropy is conserved, both normalizations are related 
by a constant. However, an inconvenient aspect of the quantity Yx is its dependence on 
the entropy degrees of freedom, gs{'t*)/9s{'to), and on possible entropy production between 
ti, and ^0- 

The final baryon asymmetry is conveniently expressed in terms of the baryon-to-photon 
ratio TjBo — nB{to)/n^{to), to be compared with the CMB measurement (2). The predicted 
value of rjBo is obtained from by accounting for the dilution factor / = N^{to) > 1, 



In the simple case of constant entropy one has / = gs/g% with g^ = 2 + 21/11 ~ 3.91. 
Assuming at the standard model degrees of freedom with a single Majorana neutrino 
in addition, one obtains g^ = 434/4, and therefore tibq — 0.036A^bo- 

In the decays of the heavy Majorana neutrinos an asymmetry in the number of lepton 
doublets, and therefore in B — L, is generated. In the following we shall sum over the three 
lepton numbers L^, and L^-. Because of the large neutrino mixings suggested by the 
solar and atmospheric neutrino anomalies, we expect this to be a good approximation. A 
refined analysis can be performed along the lines discussed in ref. [19]. A related problem 
is the role of 'spectator processes' [20] which change the naive sphaleron baryon-to-lepton 
conversion rate by a factor 0{1), since any generated asymmetry in lepton doublets is fast 
distributed among many leptonic and baryonic degrees of freedom in the plasma. In the 
following we shall ignore this uncertainty and use the naive sphaleron conversion factor 
for Nb/Nb-l, which in the standard model with one Higgs doublet is a = 28/79 ~ 0.35 
[21]. The baryon-to-photon ratio today is then given by 



(5) 



Vbo 



~ 0.013 



B-L ■ 



(6) 
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The Boltzmann equations for the time evolution of the number of heavy Majorana 
neutrinos, A^i, and oi B — L number, Nb-l, are given by [8, 9], 



^ = -{rn + Ts){N^,-N^^J (7) 

= -e^rn{N^,-K\)-rwNB-L. (8) 
Here the rate To accounts for decays and inverse decays {z — Mi/T), 

where Ki and K2 are Bessel functions, and h is the Dirac neutrino Yukawa matrix (cf. sec- 
tion 3); the inverse decay rate is given by Fjo — (rf^^/rii) Fd. 

The A^i scattering rate involves processes with the Higgs field in t- and s-channel, 

= 2r(^^) + 4r(^^^ . (lo) 

Inverse decays, AL = 1 processes (F^^j, F^^g) and AL = 2 processes (Fjv, F^^t) all con- 
tribute to the washout rate. 



= ( -Fin + 2ri'i + r^'i^ j + 2 C + 2C, (11) 



The quantities F- are thermally averaged reaction rates per particle X. They are related 
by f[^^ = ^i/rfx reaction densities 7j [8] which are obtained from the reduced cross 

sections ai{s/Mi), 

1 f°° 

^«(^) = ^7^~ / dxa(^i-j{x) ^Kx{z^) , (12) 

where and are the masses of the two particles in the initial state. Our calculations 
are based on the reduced cross sections given in ref . [22] . 

An important part of our analysis is an improved treatment of the AL = 2 processes 
which involve the heavy Majorana neutrinos A^j, i = 1 . . .3, as intermediate states. The 
reduced cross sections ctat and (T(Ar,t) read. 



1 



27r 



i<j 



(13) 



with 



Q-i CLj 



1 1 

+ + 



Di{x) D,{x) Di{x)Dj{x) 



+ 1 + - 
V X 



ttj — tti Dj{x) 



In 1 + 



X 



Oj — Qj Di[x) 



X 



+ 



x + aj X + 2a j 



In 1 + 



X 



In 1 + — 



(15) 
(16) 



(oj - aj){x + Oj + Oj) 



{2x + 3ai + aj) In ^1 + ^ 



— {2x + 3aj + Oj) In I 1 



X 



;i7) 



Here % = Mj/Mf, and 1/A(a:) = (a; — ai)/[(a; — a^)^ + ajCj] is the off-shell part of the Ni 
propagator with q = ai{h}h)%/{d>T^Y . 

The sum Jn + jN,t is conveniently separated in two parts. The first part comes from 
the resonance contribution oc x/Df, which is highly peaked around x — 1. This term is 
easily evaluated analytically in the zero- width limit, 



(18) 



For typical values of {h'^h)ii and Mi the resonance contribution dominates in the tem- 
perature range from T ~ Mi {z ~ 1) down to T ~ O.lAfi {z ~ 10). The remaining 
part is dominant at low temperatures, z > 10. For z ^ 1 the main contribution to the 
integrals 7Ar and 7^^^^ comes from the region a; ^ 1. Here the scattering amplitudes with 
Ni exchange are proportional to the light neutrino mass matrix 771^. To leading order in 
X for and d'(^N,t), one finds 



-fN{z > 1) ~ liN,t)iz > 1) 



3Mf 1 . . . 
^^tr(mtm.) . 



(19) 



For hierarchical neutrinos Iff' — tr (mj^mj,) ~ Am^^^, whereas for quasi-degenerate neu- 
trinos one has nf ~ 3m^, with mi ~ m2 — rris ^ Am^^^. 

Since all rates are expressed as functions oi z, it is convenient to replace time by 
z — Mi/T in the eqs. (7) and (8). This change of variables introduces the Hubble 
parameter since dt/dz — l/{Hz), with ~ 1.66 y^(Mf /Mpi)/2;^ where Qp is the number 
of energy degrees of freedom at t^,. Neglecting the small variation of the number of degrees 
of freedom during leptogenesis one obtains for the standard model with one right-handed 
neutrino Qp — 434/4. 



The kinetic equations for leptogenesis now read, 

^ = -{D + S) {N^, - N^^J , (20) 



= -s, D {N^^ - N'^J - W Nb-l , (21) 

where we have defined {D,S,W) — {rD,^S:^w)/{H z). In order to understand the 
dependence of the solutions on the neutrino parameters, it is crucial to note that the 
rates Fu, and also T^, except for the contribution AT^ = ~ ^^N,res + '^^N,f> 
are all proportional to {h^h)ii. The rescaled rates in eqs. (20) and (21) are therefore 
dimensionless functions of z, proportional to 

D,S.W-AW oc'^, AWoc'J^^; (22) 
the effective neutrino mass rhi [9] is given by 

Ml ^ ^ 

where mo is the Dirac neutrino mass matrix (cf. section 3). Eq. (22) implies that, as long 
as AW can be neglected, the generated lepton asymmetry is independent of Mi. 

Since AW increases with Mi, it does become important at large values of Mi. At 
z ^ 1, one easily obtains from (19), 

Ml frf 

AW{z»l)^a^^ , (24) 

z 

with = C{3)7r^ giv^H{z = 1)/Mf. 

For z < 1, AW is no longer proportional to Mim^, but depends on the heavy neutrino 
masses Mj and on the specific structure of the h matrix. An approximate upper bound 
on Fw is given by the sum of the resonance contribution F^^^^es AF{z ^ 1) taken at all 
values oi z. In this way the relativistic suppression of A^i exchange at high temperatures is 
underestimated. Analogously, a lower limit F^;^ can be obtained by treating in the off-shell 
part AF the heavy neutrinos N2 and kinematically hke A^i, which overestimates the 
relativistic suppression of A^2,3 exchange at high temperatures. The two branches F^ are 
shown in figs. la-4a. As expected, the uncertainty of Fw{z) at small z is not important 
for the final baryon asymmetry in most cases. Numerically, we find that is negligible for 
heavy neutrino masses Mi < 10^^ GeV (0.1 eV/m)^. 

Successful leptogenesis requires a departure from thermal equilibrium for the decaying 
heavy Majorana neutrinos. Furthermore, at the same time, washout processes must not 



be in thermal equilibrium. The corresponding naive out-of-equilibrium conditions are 
(Td + ^s)\z=i < H\z=i ^i^l2=i < ^L=i- These conditions are fulfilled for a typical 
set of parameters Mi = 10^° GeV, fhi — 10~^ eV, m — 0.05 eV, for which the three rates 
and the Hubble parameter are shown in fig. la. The generation of the B — L asymmetry 
for these parameters is shown in fig. lb for \ei\ = 10^^. The figure also demonstrates 
that the Yukawa interactions are strong enough to bring the heavy neutrinos into thermal 
equilibrium. The resulting asymmetry is in accord with observation. 

Increasing rhi by two orders of magnitude to rhi = 10~^ eV increases all rates while 
leaving the Hubble parameter unchanged. As fig. 2a shows, the out-of-equilibrium condi- 
tion for A^^i decays is now no longer fulfilled, and the final B — L asymmetry is reduced 
by two orders of magnitude. Increasing Mi to 10^^ GeV increases AH^ by five orders of 
magnitude, which now dominates the shape of the washout rate Vw (fig- 3a). Although 
the Yukawa interactions arc strong enough to bring the heavy neutrinos into thermal equi- 
librium, the washout rate is now so large that the final B — L asymmetry is reduced by 
three orders of magnitude. Finally, reducing fhi to 10^^ eV while keeping Mi = 10^*^ GeV 
fixed (fig. 4a), one becomes dependent on the initial conditions. The Yukawa interac- 
tions are no longer strong enough to bring the heavy neutrinos into thermal equilibrium 
(fig. 4b). Starting from zero initial abundance, iV^^ = 0, the final B — L asymmetry is 
reduced by one order of magnitude compared to fig. lb. On the other hand, assuming 
initially a thermal distribution, which may have been generated by other interactions, the 
final asymmetry is enhanced by one order of magnitude compared to fig. lb. This is the 
'way-out-of-equilibrium' case, where washout effects can be neglected, and the final B — L 
asymmetry is given by 3£i/4. 



3 The CP asymmetry 

Consider now the standard model with right-handed neutrinos. The neutrino masses are 
obtained from the lagrangian, 

jC-m = hijluJ^Rj(f) + ^MijV%-i/Rj + h.c. . (25) 

Here M is the Majorana mass matrix of the right-handed neutrinos, and the Yukawa 
couplings h yield the Dirac neutrino mass matrix itld = hv after spontaneous symmetry 
breaking, v = {(f)). We work in the mass eigenstatc basis of the right-handed neutrinos 
where M is diagonal with real and positive eigenvalues Mi < M2 < M3. The seesaw 
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mechanism [3] then yields for the hght neutrino mass matrix, 

= -niD^mJ) , (26) 

where higher order terms in 1/M have been neglected. 

The mass matrix rrii, can again be diagonalized by a unitary matrix C/^*^^ , 



/ mi \ 

m2 
V ms J 



= -Dm , (27) 



with real and positive eigenvalues satisfying mi < m2 < rris Inserting eq. (26) one finds, 

v^U^^^^hD^^^h^U^'^ = Dm , (28) 

which means that 

n = vD-'/'U^'^^^hDl,'/' (29) 

is an orthogonal matrix, flff" — I [23]. This also implies lm(0^0)ii = 0, from which one 
immediately obtains 

-llm(C/MtA)^, = -5:^Im(C/«''>tA)^ (30) 

'"'1 . '"'i 

The CP asymmetries in the decays of the heavy Majorana neutrinos arise at one- loop 
order from the interference of the tree level amplitude with vertex and self-energy correc- 
tions [24, 25, 26]. In the following, we shall restrict ourselves to the case of hierarchical 
heavy neutrinos, i.e. Mi ^ 1^2,^3. The CP asymmetry ei for the decay of A^^i is then 
easily obtained by first integrating out the heavier fields N2 and N^. From the tree am- 
plitude and the effective lepton-Higgs interaction one then obtains the useful expression 
[18], 

where corrections 0{Mi/M2,3) have been neglected. Using eqs. (27) and (30) one easily 
derives, 

_ 3 Ml ^ Am?. 1°' (''") 
where Am^^ = m? — mf, and 

h = U^'^^h (33) 
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is the matrix of Yukawa couplings in the mass eigenstate basis of hght and heavy Majorana 
neutrinos. In the interesting case M2 — Mi — 0{Mi) the CP asymmetry ei is enhanced 
by a factor 0{1). However, one then has to study the decays of both heavy neutrinos, A^i 
and N2 [16]. 

In the case of hierarchical neutrinos, ~ (Am^^^)^/^ ^ m2 — (Am^^J^/^ ^ mi, 
eq. (32) yields an upper bound on the CP asymmetry [27] 



3 MijAmlJ'/' 



where we have neglected the term oc Am^^;. For an inverted hierarchy ~ m2 — 
(Am^^^)^/^ 3> mi and ml — m| = Am^gj both terms contributing to eq. (32) are approx- 
imately equal to (Am^j^)^/^, and the upper bound is larger by a factor two. Further, 
for Yukawa couplings with \hij\ < {h^^l = C(l), one has M3 ~ v'^/m^. The maximal CP 
asymmetry is then a measure of the hierarchy among the heavy Majorana neutrinos, 

|£i|"^^~0.1^ . (35) 
' " M3 ^ ^ 

It is remarkable that the upper bound (34) is frequently saturated in models with hi- 
erarchical neutrino masses (see, e.g., [19], [28]- [32]). One reason is that in the leptonic 
mixing matrix U = q\\ elements, except f/gS) ai'c known to be 0{1). This is 

often explained by the structure of U'^'^\ whereas U^^^ has small off-diagonal elements. 
The Yukawa matrix = U^'^^'^hi, has then naturally large off-diagonal elements, even if 
h,y is almost diagonal. For hierarchical neutrinos the upper bound on the CP asymmetry 
is then easily saturated. 

In terms of the quadratic mean m/ \/3, with 



m = \j rn\ + ml + m| , (36) 
the bound (34) may be written as 



I I . 3 Mim 

ki < — — ■ 37 



For quasi-degenerate neutrinos, where mi ~ m2 — m^ ~ m/\/3 3> (Am^^^)^/^, one gets 
the stronger bound 

< bUB^ . (38) 



There is no well defined boundary between hierarchical and quasi-degenerate neutrinos. In 
our analysis we shall choose m — l eV, such that fn/ -\/3 ~ 0.58 eV ^ •>/ Am^^^ ~ 0.05 eV. 
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The interesting case where all neutrinos have masses nii — C(O.leV) will be discussed 
elsewhere [16]. 

What is the allowed range of the masses fhi and Mi which are crucial for leptogenesis? 
Assuming Yukawa couplings \hij\ < 1, one has 

Mi<M3^ — . (39) 
ma 

The effective neutrino mass fhi is given by (cf. (29)), 

^i = ^EI^'iHE^^I^'iI- (40) 

i i 

From the orthogonality of fl one obtains the lower bound [33] , 

^1 > ^ > mi ^Re(Q^i) = mi . (41) 

i i 

The orthogonality condition for fl reads explicitely, 

Re(Jl^Q)ii = ^Y1 = 1 ■ (42) 

Unless there are strong cancellations due to phase relations between different matrix 
elements, one then obtains 

^ "^3^ E - "^3^ E = ma . (43) 

i i 

In our analysis we shall therefore emphasize the range mi < mi < m^. 



4 Constraints on neutrino parameters 

In this section we shall compare the predicted baryon asymmetry with the value mea- 
sured from CMB (cf. (2)), using the relation (6). Incorporating also the bound on the CP 
asymmetry described in the previous section will enable us to derive an allowed region 
for the three parameters Mi, mi and, very interestingly, m. Before describing the numer- 
ical results it is useful to discuss some general properties of the solutions of the kinetic 
equations. 

The abundance N^i (z) of the heavy Majorana neutrinos is independent of the asym- 
metry Nb-l- Thus, for a given solution Njsi^^z), or equivalently A (2;) = Nj^^{z) — N'^^{z), 
the solution for the asymmetry Nb-l{z) reads, 

Nb-l{^) = ArB_L(^in)e-^^n -^^'^(^'^ + -£i«(^;mi,Mi,m) , (44) 
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where 2;in <^ 1, and the efficiency factor k is given by 

^[z) = \ r dz' D{z')A{z') e-/;"^^"^(-") . (45) 

Since we assume that the initial asymmetry is zero, Nb-l{z) is proportional to the CP 
asymmetry ei. 

A interesting, limiting case is 1^ = and — 0, i.e. A = A^jvi, which corresponds 
to far out-of-equilibrium decays with vanishing washout. The final baryon asymmetry 
then takes the form 

7VO_^ = £i7V^,(^,), (46) 

where za is some temperature at which the heavy neutrinos are out of equilibrium and 
have not yet decayed. For an initial thermal abundance one then obtains A^B_i = 3£i/4, 
and therefore kq — 1. Within thermal leptogenesis this is the maximal asymmetry, and 
one always has Kq < 1. 

In special parameter regimes it is possible to obtain explicit analytical expressions 
for the efficiency factor ^(z; mi. Mi, m) [19]. For small fhi and fn washout effects are 
small. Starting from A^^'^ = 0, the final asymmetry is then proportional to the generated 
Ni abundance. From the reaction densities given in [22], one easily obtains at 2; <^ 1: 
D oc Mpifhiz^ / v'^ and S oc Mpifhi/v'^. Hence, one obtains from eq. (20), 

K N^.iz,) oc . (47) 

Also interesting is the case of large fhi and small Mi, where the washout is dominated 
by the W — AW term. At large temperatures, z < 1, the thermal A^i abundance is 
then quickly reached. At small temperatures, z > 1, the Ni abundance is reduced by the 
AL = 1 scatterings and by the decays. The reaction densities '~fs,D decrease exponentially, 
and the efficiency factor is determined by the A''i abundance at freeze-out z [19], 

K oc (z) oc oc . (48) 

ls,D mi 

For intermediate temperatures, i.e. z ^ 1, the Boltzmann equations have to be solved 
numerically. 



4.1 Numerical results 

In fig. 5 we have plotted^ the predicted present baryon asymmetry as a function of the 
parameter fhi for different values of Mi, assuming a typical value of the CP asymmetry, 

^AU numerical results have been crosschecked by two independent codes. Figs. 1-4 and 5b-7b are 
based on one code; the other code has been used for figs. 5a- 7a. 
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£i = — 10~^. We have calculated the baryon asymmetry for two values of m. In fig. 5a, 
fn — 0.05 eV ~ ^/Am^^, which corresponds to the case of hierarchical neutrinos; m — 
1 eV in fig. 5b, which represents the case of quasi-degenerate neutrinos with all hght 
neutrino masses approximately equal, mi = l/\/3eV ~ 0.58 eV. We have performed the 
calculations both for an initially vanishing Ni abundance (thick lines) and for a thermal 
initial A^i abundance, iVj^^ =3/4 (thin lines). 

Following [19], we show in fig. 6 the dependence of the baryon asymmetry on the 
parameters rhi and Mi by means of Iso-Kq curves of the efficiency parameter. Figs. 6a 
and 6b correspond to the cases m = 0.05 eV and m = 1 eV, respectively. For vanishing 
initial A^i abundance the enclosed domains have a finite extension in rhi] for thermal 
initial abundance there is no boundary at small rhi. 

Figures 5 and 6 clearly show the existence of two difi'erent regimes: the domain of 
'small' iVi masses. Mi < 10"*^^ GeV(0.1 eV/m)^, and the domain of 'large' A^i masses. 
Ml > 10"*^'^ GeV(0.1 eV/m)^. The dependence of the boundary on fn is determined by 
the behaviour of the non-resonant washout rate, AW oc Mifn'^. Note, that in previous 
studies of the washout effects [10] the obtained Mi dependence was a result of an assumed 
behaviour of Al^ at z ^ 1. As our discussion in section 2 shows, such an assumption is 
unnecessary, since the behaviour of AW is governed by Mifff. 

4.1.1 Small Ml regime 

In this case the non-resonant part AW of the washout rate is negligible in first approxima- 
tion, and there is only a small dependence on Mi and fh in the final asymmetry. Hence, 
the efficiency factor Kq depends approximately only on the parameter rhi. From fig. 5 
one reads off that in the case of zero initial A^^i abundance the final asymmetry reaches 
its maximum at fh\^'^^ ~ 6 x 10^"^ eV. We can then express the baryon asymmetry in the 
following form, 

VBo = l.5xlO 9 (49) 

where Kq^^^ = Ko{rhl'^^^) ~ 0.16 is the maximal efficiency factor kq obtainable for zero 
initial A^i abundance. 

The behaviour of the efficiency factor for small and large values of rhi is known analyt- 
ically (cf. (47), (48)). Approximately, one has oc fhi for fhi <^ m^^^ , and kq oc 1/fhi 
for fhi ^ fhi'^^^. A good fit for «;o(^i) is given by the following expression, 

Ko{fhi) ~ 0.24 (x_ e-''- + x+e-''+) , (50) 
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with 

^) . (51) 

m±J 

and m_ = 3.5 x 10^"^ eV, m+ = 8.3 x 10^'^ eV, a = 0.1. The corresponding baryon 
asymmetry is shown in fig. 5a; the fit is optimal for Mi = 10^ GeV. 

In case of thermal initial Ni abundance the maximal efficiency factor is obtained in 
the unphysical limit rhi 0, where the heavy neutrinos decouple completely from the 
thermal bath. One then has 

r^^o^" ~ 0.96 X lO^Vi . (52) 

On the other hand, for sufficiently strong coupling of the heavy neutrinos to the thermal 
bath, i.e. rhi > rhi'^^^, there is no dependence on the initial neutrino abundance, and 
the two cases N^^ = and iV]^^ =3/4 give the same final baryon asymmetry. In this 
case possible asymmetries generated in the decays of the heavier neutrinos N2 and A^3 
have also no effect on the ffnal baryon asymmetry. Note, that in comparison with the 
analytical results of ref. [19], we obtain an efficiency factor k which is about three times 
smaller. 

4.1.2 Large Mi regime 

In this regime the washout rate AW oc Mi dominates and the Iso-kq curves in the 
(m-i, Mi)-plane tend to become independent of rhi with increasing Mi. However, some 
dependence on rhi remains due to the effect of D and S (cf. (20), (21)) on the Ni abun- 
dance. This is why the curves never become exactly horizontal. Note, that the curves 
tend to converge in such a way that a small variation of Mi leads to a large variation of 
the final asymmetry. In this region any small change in the kinetic equations can have 
large effects on the final asymmetry. On the other hand, a curve corresponding to a given 
value of the final asymmetry, e.g. tJbo^^, is rather insensitive to small changes of the 
kinetic equations. Thus in this region the uncertainty of AW at small values of z is, 
fortunately, not important. We therefore conclude that, within our framework of thermal 
leptogenesis, we have obtained a description of the baryon asymmetry in terms of just 
four neutrino parameters: Si, Mi, rhi and fh. 

4.2 CMB constraint 

The CMB constraint on neutrino parameters is given by the requirement 

r)Bo ^ 0.96 X 10-2 \ei\ Ko{mi, Mi,m) = 77^^^ . (53) 
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In the case N^^^ — 0, the substitution of the efficiency factor by its maximum Kq"^^^ yields 
an important lower limit on the CP asymmetry £i, 

hi > 4.0x10-^ f lO-io ) ^ 2.4 xlO-^ (54) 

Here the last inequahty corresponds to the {3a) lower limit for rj^^^ (cf. (2)). For small 
Ml, the CMB bound becomes a determination of fhi as function of the CP asymmetry. 
There are obviously two solutions: m]" < 'm\^^^ and mf > m^^^^. Prom the fit (50) one 
easily finds fnf for CP asymmetries sufficiently far above the bound (54), 

/ CMB ^r)-6\l-l 

r - ^•Sxl''-=ev(^^) , (55) 



m 



6 X 10-^0 |£i| 



0.9 



mt ^ 2.8 X 10-3 eV ^^l^^g^^j . (56) 

As Si decreases to the lower bound (54), m]~ and approach m^'^^^ ~ 6 x 10^^ eV. 

In the case of a thermal initial A^^i abundance the bound (54) gets relaxed by a factor 
l/A;^'^ ~ 6.4, and one obtains the usual bound for kq = I in eq. (53), 

kil > 6.3 X 10-« (^^^) > 3.8 X 10-« . (57) 

This bound can only be reached for mi <S fhi'^'^^, i.e. very small values which are not 
easily obtained in models of neutrino masses. Note, that rhi > mi, the smallest neutrino 
mass eigenvalue [33] and that rhi^^^ ~ O.ly^Am^. Purther, one has to worry about the 
production mechanism of a large initial abundance of extremely weakly coupled heavy 
neutrinos. 



4.3 CMB constraint plus CP bound 

Until now we have treated the CP asymmetry Si as an independent parameter. However, 
as discussed in section 3, for hierarchical as well as quasi-degenerate neutrinos, the CP 
asymmetry satisfies an upper bound |£i| < e{Mi,m) (cf. eqs. (37), (38)). Together with 
the CMB constraint (53) this yields the following restriction on the space of parameters 
rhi, Ml, and m, 

r/go"(mi,Mi,m) ~ 0.96 x IQ-^ e{Mi,m) Ko{mi,Mi,m) > rg^"" . (58) 
Por hierarchical neutrinos eq. (37) yields the upper bound on the CP asymmetry, 

16 



Combining this with the CMB constraint (54) yields a lower bound on Mi [27] . For zero 
(thermal) initial A^i abundance one obtains, 



In the case of inverted hierarchy the upper bound on the CP asymmetry is twice as large 
and therefore the lower bound on Mi is twice as small. 

In the case of quasi-degenerate neutrinos eq. (38) implies a stronger bound on the CP 
asymmetry, 



The corresponding lower bound on Mi reads for zero (thermal) initial A^i abundance. 



The bounds (59) and (61) on the CP asymmetry seem to suggest that by increasing 
Ml any value |£i| < 1 can be reached. One may therefore expect it to be rather easy to 
satisfy the CMB constraint (53). However, in the regime of large Mi the efficiency factor 
/«o(toi. Ml, m) is exponentially suppressed. This dominates the linear increase of |£i|™^ 
with Ml and leads to an upper limit on Mi. 

This situation is illustrated in figs. 7a and 7b, where \so-rf^Q^ curves in the (mi, Mi)- 
plane are shown for m = 0.05 eV and fn = 1 eV, respectively. In the first case we have 
assumed a normal hierarchy. In the second case we have used |£i| < lO^^Mi/lO^'^GeV 
corresponding to Am^^^ ~ 2.9 x 10~^eV^. For values rhi and Mi enclosed by these curves 
the baryon asymmetry tj^q^^ can be obtained. In the case of N^^_^ — 0, the allowed regions 
are closed domains; for N^^_^ =3/4 there is no lower limit on rhi. 

As discussed in section 3, the effective neutrino mass rhi is bounded from below by 
mi, the mass of the lightest neutrino, and also likely to be smaller than m^, the mass 
of the heaviest neutrino. For quasi-degenerate neutrinos this implies mi ^ m/ ^/S. This 
value is indicated by the vertical line in fig. 7b from which one reaches the conclusion 
that quasi-degenerate neutrinos are incompatible with leptogenesis. A similar conclusion 
was reached in ref. [33] based on results of [10], assuming zero initial A^i abundance. 
As our analysis shows quasi-degenerate neutrinos are strongly disfavoured by thermal 
leptogenesis for all possible CP asymmetries and initial conditions. Possible ways to 
evade this conclusion are a resonant enhancement of the CP asymmetry in the case of 




(60) 




(61) 




(62) 
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heavy neutrino mass differences of order the decay widths, |M2,3 — Mi\ — 0(rj) [34], or 
a completely non-thermal leptogenesis [33] . 

In the case of hierarchical neutrinos there is a lower bound on Mi as function of mi, 
which can be read off from fig. 7a in the case of a normal hierarchy, while for inverted 
hierarchy this bound is twice as small. No upper bound stronger than M3 ~ 10^^ GeV 
exists. For thermal initial conditions all values of fhi < 777,3 ^-re allowed. 

Comparing figs. 7a and 7b, it is evident that a more stringent upper limit than 1 eV 
exists for fn. The determination of this precise bound goes beyond the goal of this paper 
and is left for future work [16]. 

Any point in the space {fiBOjmi, Mi) which lies below the surface ^7^0^ and above the 
plane tJbo^^ represents a possible set of neutrino parameters with some value < 
As discussed in section 3, in many models one has \ei\ ~ The CMB bound then 

yields a precise relation between rhi and Mi. In particular, for the interesting range 
O.leV > fhi > rh^^^^ ~ 6 x 10~^eV, one has the simple relation. 

This corresponds to the scenario originally proposed in [28], with rhi ~ -y/Am^ ~ 
0.005eV, where we have assumed the LMA solution [35]. The corresponding baryogenesis 
temperature is ~ Mi ~ 10^° GeV. 



5 Conclusions 

We have studied in some detail the minimal version of thermal leptogenesis. In this frame- 
work the decays of A^i, the hghtest of the heavy Majorana neutrinos, are the source of 
the baryon asymmetry, and only Yukawa interactions generate the initial heavy neutrino 
abundance; further, the connection between light and heavy neutrinos is given by the see- 
saw mechanism. The final baryon asymmetry is then determined by only four parameters: 
the CP asymmetry £1, the heavy neutrino mass Mi, the effective light neutrino mass rhi 
and the quadratic mean m of the light neutrino masses. 

The constraint from the cosmic microwave background on the baryon asymmetry 
strongly restricts the allowed range of neutrino parameters. For small values of Mi, 
the efficiency factor k depends only on fhi. Together with the upper bound on the CP 
asymmetry, this yields an important lower bound on Mi and therefore on the baryogenesis 
temperature, Tg ~ Mi > 2.4 x 10*^ GeV. 
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For large values of Mi, washout effects, proportional to Mifrf, become important. 
This leads to the conclusion that the case of quasi-degenerate neutrinos, fn » y^Am^, 
is incompatible with thermal leptogenesis. On the other hand, in the case of hierarchical 
neutrinos, the full range of light and heavy neutrino masses considered in grand unified 
theories is viable. The intermediate regime, = (9(0. leV), remains to be studied in 
detail. 

Global fits of solar and atmospheric neutrino data indicate large mixing angles for 
neutrino oscillations in both cases. With respect to neutrino masses, leptogenesis favours 
the possible case of a mild hierarchy, 1712 — a/Attt^ ~ 0.1 \/Am^^i ~ 0.1 ms. Due 
to the large mixing angles, the a priori possibility mi <^ 1712 appears unlikely. With 
mi < fhi < m3, one then arrives at the conclusion mi ~ 0.1\/Am|^. Furthermore, 
in many of the recently discussed models for neutrino masses the CP asymmetry Si is 
close to its upper bound. This finally determines the heavy neutrino mass Mi and the 
baryogencsis temperature to be Tb ~ Mi = O{10^^ GeV). 

The unfortunate prediction of thermal leptogenesis is that the search for absolute 
neutrino masses in tritium /5-decay and CMB combined with Large Scale Structure will 
be unsuccessful in the near future. Our present bound, ^^mj < \/3 eV, can still be 
significantly improved [16]. However, the search for neutrinoless double /3-decay with a 
sensitivity |mee| = 0(v^Am|^) [36] may very well be successful. 
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figure 1a 
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Figure 1. Mi = lO^'^ GeV, mi = 10~^ eV. (a) Rates normalized to the expansion rate 
at z = 1. The two branches for Tw at small z correspond to the upper (lower) bounds 
i^w) (^^® text), (b) Evolution of the A'^i abundance and the B — L asymmetry for 
El = — 10~^ and m = 0.05 eV, both zero and thermal initial A'^i abundance. 
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figure 2a 
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figure 2b 
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Figure 2. (a) Rates and (b) evolution of A^i abundance B — L asymmetry; parameters 
in figure 1, except mi = 10~^ eV. 
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figure 3a 
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figure 3b 
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Figure 3. (a) Rates and (b) evolution of A^i abundance B — L asymmetry; parameters 
in figure 1, except Mi = 10^^ GeV. 
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figure 4a 
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figure 4b 
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Figure 4. (a) Rates and (b) evolution of A^i abundance B — L asymmetry; parameters 
in figure 1, except rhi = 10~^ eV. 
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Figure 5. Predicted baryon-to-plioton ratio rjBo as a function of mi, for ei = — 10~ 
and for the indicated values of Mi. The horizontal sohd and dashed hues indicate the 
mean and the upper/lower values (3cr) of 77^9^^, respectively (see eq. (2)). (a) Hierarchical 
neutrino case (m = 0.05 eV ~ a/ Am^^j^). The squares denote the fit (50). (b) Quasi- 
degenerate neutrino case {m^. ~ 0.58 eV). 
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figure 6a 




figure 6b 
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Figure 6. Iso-kq curves, (a) Hierarchical neutrino case (m = 0.05 eV ~ 
Quasi-degenerate neutrino case {m^. ~ 0.58 eV). 
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figure 7a 
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figure 7. Iso-?7^Q^ curves. The region within the dashed hnes is currently allowed by 
the CMB constraint (cf. 58) for the {3a) lower rj'^Q^^ value, (a) Hierarchical neutrino case 
(m = 0.05 eV ~ a/AtttI^). (b) Quasi-degenerate neutrino case (m^. ^ 0.58 eV). The 
dashed vertical line and the arrows refer to the bound (41). 
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Erratum 



As pointed out by Giudice et al. [1], a proper treatment of the real intermediate state 
contribution to the AL = 2 scattering processes also subtracts the term given in Eq. (18). 
This corrects a result arising from the subtraction procedure introduced in [2], which 
has been extensively used in the literature. Removing the spurious term (18) reduces 
the washout rate T\y by a factor of about 2/3. As a consequence, the numerical factors 
change in Eqs. (60) and (63): 



A detailed discussion can be found in [3]. 
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